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O ' Abstract 

We write the thermodynamic Bethe ansatz for the massive OS'p(2|2) Gross Neveu and 
, sigma models. We find evidence that the GN S matrix proposed by Bassi and Leclair [12] 

^ ' is the correct one. We determine features of the sigma model S matrix, which seem highly 

5^ , unconventional; we conjecture in particular a relation between this sigma model and the 

complex sine-Gordon model at a particular value of the coupling. We uncover an intriguing 
duality between the OS'p(2|2) GN (resp. sigma) model on the one hand, and the S'0(4) sigma 
(resp. GN model) on the other, somewhat generalizing to the massive case recent results on 
OS'p(4|2). Finally, we write the TEA for the (SUSY version of the) flow into the random 
bond Ising model proposed by Gabra et al. [39j , and conclude that their S matrix cannot be 
correct. 



1 Introduction 

The study of 1 + 1 quantum field theories with supergroup symmetry is a lively and difficult 
topic, central to several key areas of modern physics, in particular the solution of critical points 
in non interacting disordered systems, and the AdS/CFT conjecture. 
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It is natural to try to make progress by focussing on integrable theories. By analogy with 
the ordinary case, one can expect the simplest possible situations to be encountered [IJ in the 
0Sp{n\2m) Gross Neveu and sigma models. Yet, apart from the special case of OSp{l\2) which 
behaves much like an ordinary group [213], only little progress has been accomplished, largely 
for what seems to be technical reasons. Another exception is the case of OSp{4:\2), where, among 
other fascinating properties [3], a striking duality between the GN and sigma model was uncovered 

mm- 

The case of OSp{2\2) is particularly tentalizing, since it is the 'simplest' after 05^(112), and 
plays a major role in disordered systems. Yet, even here, significant progress has only occurred 
up to now in the critical case, where WZW theories have been solved [7] , and their relation with 
spin chains understood ^ ^ . Our main goal in this paper is to understand the massive theories: 
the GN and OSp{2/2)/OSp{l/2) 'supersphere' sigma models. 

To proceed, we recall some general features in the ordinary 0{N) bosonic case. The actions 
are well known [lOJ. The S matrices involve generically particles in the (vector) representation. 
Introducing the useful graphical representation of invariant tensors there are generically two 





Figure 1: Graphical representation of the invariant tensors appearing in the 5* matrix. 

known models whose scattering matrix for the vector representation has the form in figure [TJ 
with none of the (Tj's vanishing. They are given by 

_ 2i7r 

~ ~{N -2){i^-e) 

2i7r 

0-3 = -TTT ^ CF2 (1) 



{N -2)9 



with two possible choices for a2'- 



N-2 ' 2iTT V 2 N-2 2iTT 



The factor ex J does not have poles in the physical strip for > 0, and the corresponding S matrix 
for A > 3 is believed to describe the 0{N) /0{N — 1) sphere {S^~^) sigma model. The factor 
(72 does not have poles in the physical strip for A < 4. For A > 4, it describes the scattering 
of vector particles in 0(A) Gross Neveu model. Recall that for A = 3,4 the vector particles in 
the GN model are unstable and disappear from the spectrum, that contains only kinks. Some of 
these features are illustrated for convenience in figure 2. 

Note that at vanishing rapidity, the scattering matrix reduces to S{9 = 0) = ^I. This is in 
agreement with the fundamental particles being bosons in the sigma model , and fermions in the 
Gross Neveu model [H]. 

It was argued in [1] that for A negative and |A| large enough, a similar S matrix - where 
now the invariant tensors have to be interpreted within 0Sp{m\2n) representation theory - still 
describes the GN and sigma models provided one sets A = m — 2n. The detailed calculations 
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region where a"*" has no poles 



region where a has no poles 

Figure 2: Pole structure of a2 as a function of A^. Except for a finite region at small A^, the 
generic S matrix with ctj^ describes the sigma model and the one with the GN model. 

performed in \Jl show that one the roles of the two models remains the same, that is, cr+ still 
corresponds to the sigma model and (T_ to the GN model. 

This has several interesting consequences. One is, that the GN model exhibits bound states 
while the sigma model does not, for N large, ie m >> 2n; meanwhile, the sigma model does 
then exhibit bound states - and the GN none - for — large, ie m << 2n. Where exactly 
these behaviours may be modified at small is a partly open question. We emphasize that, 
although the S matrices are presumably known in the generic case, no systematic check has been 
performed, as writing the TBA has proven, up to now, too cumbersome. 

A crucial observation based again on [1] and also [3] is that the massive continuum limit of the 
integrable lattice model based on the fundamental of OSp{m\2n) always exhibits bound states. 
Therefore it necessarily corresponds to the GN model for N > 2 and the sigma model for N < 2 
(including A^ negative of course). The case N = 2 being at the border suggests that the sigma 
and the GN model are then dual of each of other. This is of course well known for 0(2), and 
has recently been established, for a certain choice of boundary conditions, in the orthosymplectic 
case [5], [6]. Note that the S matrices for the sigma and GN model at A^ = 2 also coincide. They 
describe a special point on the critical lines of either (dual) model. 

The case A^ = {OSp{2n\2n)) is the next case that seems really worth studying. Indeed, right 
at A'^ = 0, (T^ = a2 (and the foregoing S matrix coincides with the one in |12j). while it is clear 
that the GN and sigma models are different (see below for detailed actions etc) . Something is thus 
missing from the general picture. On the other hand, the TBA should not be too complicated. 
It is particularly tractable for a special choice of Dynkin diagram first studied in |13j. 

The OSp{2\2) has also the potential to exhibit interesting duality properties. To see this, 
consider another - a priori more anecdotical - remark about the beta functions of the sigma and 
GN models 0. In an expansion in the coupling constant the first two non trivial coefficients obey: 

—77 = — Sigma model 

/32 A^ - 2 

| = -^GNmodel (3) 

These two results are exchanged under N — 2 — )■ 2 — N. This suggests there might be some 
relationship between the 0{N) GN model and 0(4 — A^) sigma models, where an OSp supergroup 
has to be used whenever necessary. 

Of course, these two models cannot be related in general: among other things, their UV 
central charges don't match for arbitrary N. Nevertheless, as observed in some of the S 



^Note that the /3 functions do not exhibit any singularity when going from positive to negative values of A'^. 
The beta function of say the GN model for A'^ negative is simply obtained from the one for the usual GN model by 
substituting a negative value of A'^ [14]. Similarly, the /3 functions are not affected by the "twisting" of the fermion 
boundary conditions. 
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matrix elements match exactly under this transformation. One has indeed, 

[at - 4] (N) = [a^ + ^3^] (4 - N) (4) 

This implies that, at zero temperature and with the proper applied field h , the ground state 
energies of the 0{N) GN model and the one of the 0(4 — N) sigma model do indeed coincide 
as functions of h/M (and in particular, their zero field energies - the regularized version of the 
ground state energy as a function of the mass M - coincide). 

Now, if we take the 05^(212) GN (A^ = 0) model we have c = 0,c^q = 31, while for the 
50(4) sigma model we have c = 3 (no fermions in this case). So some identificatiom might be 
possible. Similarly, if we take the 05*^(212) sigma model we have c = —1, Cgg = 2 while for the 
50(4) GN model, which is a sum of two independent sine Gordon models at the marginal point, 
c = 2, again allowing a potential identification. This pattern generalizes immediately to the case 
of OSp{2n\2n) and OSp{2n\2n - 4): 

OSp{2n\2n) GN model : c = 0, Cgg = 3n 
0Sp{2n + 2|2n — 2)Sigma model : c = 3, c^q = 3n 
05^(271, 1 2n) Sigma model : c = —1, Cg|f = 3n — 1 
OSp{2n + 2\2n - 2)GN model : c = 2, c^ff = 3n - 1 (5) 

One of the questions we will investigate in this paper is the potential duality between these two 
families of models for the simplest case n = 1. 

Our strategy will be based on non perturbative calculations of the free energies at finite 
temperature of all models using an S matrix and thermodynamic Bethe ansatz approach. 

Because of the common underlying osp{2\2) = sl{2\l) symmetry, the different Bethe ansatze 
involved have a common structure of roots and basic equations, which coincide with those for 
the integrable spin chains. We thus discuss these common features first. 

2 Solving the monodromy problem: the Bethe ansatz for inte- 
grable SL(2|1) vertex models and spin chains 

2.1 The Bethe equations in the fermionic grading 

We start with the Bethe equations for the sZ(2/l) spin chain built out of alternating 3 and 

3 representations, that is the Hilbert space is 3®^ (8> 3*^^. There are different Bethe ansatze 
available, and as often, chosing the right one can be a big help. We take the Bethe ansatz based 
on a choice of purely fermionic simple roots |l5j for the s^(2|l) algebra. 

(8«8) 



Figure 3: Dynkin diagram leading to our Bethe equations 

^We denote in general by Cgg- the central charge in the (Neveu Schwarz) sector where the fermions have an- 
tiperiodic boundary conditions (on the cylinder and on the plane, since their dimension is integer. 
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Recall that the Bethe equations read then [131 E] 

n,+^Y ^^u^-j^^ , = l,...,iV 
Uj - 7/3 - « 

^V = n ^""^'- + \ a = l,...,M (6) 

with the quantum numbers S = ^^^^^^and = L — ^^2^ ■ The corresponding energy, for the 
regime of interest here, is is 

E oc 



A crucial observation [8] which allows one to draw some conclusions quickly is that this system 
of equations admits a particular solution made of sets such that N = M and uj = 7^. The roots 
then have to satisfy 

a system of equations that is identical with the Bethe equations for the spin 1 antiferromagnetic 
XXX chain (we often refer to the latter simply as XXX, or the Takhtajan-Babujian [16], chain). 
Notice however that the dynamics is such that the gaps in our system will be twice the gaps in 
the latter chain: sound velocities have to be adjusted appropriately. Often, we will see that our 
results are in some sense a 'doubling' of those for this XXX chain. 

The argument extends to a chain with alternating fully symmetric and antisymmetric repre- 
sentations, represented by Young diagrams with p boxes . The bare Bethe ansatz equations then 
look exactly as ([6]) but in the left hand side, the factors i are replaced by pi, and similarly for the 
energy ([7|) the 1 in the denominator is replaced by p'^. This time, the symmetric sector = M, 
Uj = Ji3, is related with the XXX chain of (integer) spin p. 

The argument also works for the model in the fundamental representation of osp{2/2). In 
that case indeed, the bare Bethe ansatz equations and the energy still read as ()6|7p but with the 
factor i replaced by i/2 (ie, formally, p = 1/2 in the foregoing discussion): 

Uj + t/2\ _ -TT lij - 7/3 + ^ 



V2y ^}uj-jp 



N 



Ja + t/2\ ^ Y]j^^_Uk + i^ a = l,...,M, fundamental of 05p(2 1 2) (9) 

This time the symmetric sector is related with the spin 1/2 chain. Note that, while the equations 
([6]) hold for a system made of L representations 3 and L representations 3, when it comes to the 
fundamental of osp{2/2), Q holds for L representations 4. 



2.2 The solutions for large chains 

Like for the XXX chain, the types of solutions of the Bethe equations do not depend on the source 
terms, and can be discussed in full generality H. The standard way to classify the solutions is 

■^The following is borrowed from ^ with only small modifications. 
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to consider configurations of spectral parameters that sit on poles of the bare scattering kernels 
(the right-hand-sides in (|6l9p ). A simple calculation yields the following types of "strings" 

(1) "Reals": 

unpaired, purely real spectral parameters Uj and 7/3. 

(2) "Wide strings": 

"type-I": composites containing n — 1 7's and n u's (n > 1) 



7, 



(n,n— 1) 
''a,k 
(n.n— 1) 



+ ^ + 1 - 2A;) , k = l...n 



(n,n-l) 



+ i(n-2j), j = l...n-l, 'u^''""^) G IR 



(10) 



"type-II": composites containing n 7's and n — 1 u's {n > 1) 



(n,n— 1) 
(n,n— 1) 



ui"'""^) + i (n + 1 - 2A;) , k = l...n 



u 



(n,n— 1) 



+ i(n-2j), j = l...n-l, u^^'^'-^^gU 



ill) 



Figure 4: Wide strings of lengths three and two respectively. The circles/triangles denote the 
positions of the ti's/7's involved in the string. 

(3) "Strange strings" : 

composites containing n 7's and n u's (n > 1) 



u 



{n,+) 
a,k 



or 



u 



(n,-) 



("-) 

7(3. fc 



+i(n+l-2k 



,(«) 



i(n + 1 - 2k 



u'^^^ -i{n + l-2k 



u^^^ +i{n + l-2k 



n 



A; — 1, , 

Ua G m, 

A; = 1, . . . , n , 



fl2) 



(13) 



6 



It is fundamental to observe that solutions of this type are fundamentally different from 
the usual string solutions in that the set of roots on one level of the Bethe equations is not 
invariant under complex conjugation. The other solutions discussed above (reals and wide 
strings) are invariant under this operation. This is similar to what was found recently for the 
anisotropic sl{3) chain related with the complex SU{3) Toda theory in the continuum limit 
|17| . Of course this non invariance reflects the non hermitian nature of the superalgebra 
hamiltonian: it casts doubts on analyzing the 'thermodynamics' of the system, for instance, 
but we shall see that naive calculations seem to give the correct results anyway. Note that 
although "strange strings" are not invariant under complex conjugation, the corresponding 
energy ©is still real because it depends symmetrically on the u and 7 parameters. 



II 



II 



■I 



+ 



▲ 



▲ 



Figure 5: The two types of n = 1 strange strings 



Figure 6: A strange string with n = 2. 



(4) "Narrow strings": 



composites containing n 7's and n u's (n > 1) 



7^"^ = n(r'") + ^(n + l-2j) 



j = l...n 



k = 1 . . .n 




(14) 
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Narrow strings may be thought of as special cases of strange strings or wide strings in the 
fohowing sense. 

— Combining a "+"-strange string of length n and centre u^"^ with a a "-"-strange string 
of length n and centre u^"^ we obtain a narrow string of even length 2n. This is shown 
for n = 1 in FiglTI^a). 

— Combining a type-1 wide string of length n and centre with a type-2 wide 
string of length n and centre ^("'"-"1) we obtain a narrow string of length 2n — 1. This 
is shown for the case n = 2 in Fig. EKb). 

A Im 




Figure 7: Combining a pair of 
length 2; 



'-|-" and "-" strange strings of length 1 gives a narrow string of 



Im 



Re 



Figure 8: Combining a type-1 wide string on length 2 with a type-2 wide string of length 2 gives 
a narrow string of length 3; 

It is clear from our discussion that narrow strings are not "fundamental" but are merely 
degenerate cases of strange string solutions of the Bethe ansatz equations. 
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3 TBA for the 3,3 and the 4 (fundamental of OSp(2|2)) spin 
chains 

We now derive the TBA equations for the gapless spin chains to check (some aspects of) the 
completeness of our classification of solutions, and prepare subsequent discussion. 

To start we derive equations for the densities of these various solutions in the limit of large 
systems. This is done as usual by considering that the various roots in the complexes are slightly 
off, and eliminating the corresponding small shifts through the Bethe equations, resulting in 
'higher' Bethe equations. After taking the logarithm and naively differentiating, we get a system 
which we write symbolically as 

pi + = Si + (15) 
j 

(note we have included a 2ii into the definition of the source and kernels). We introduce the 
fundamental objects 

and the Fourier transformation 

/(x) = j due'-^fiu), f{u) = i- y dxe-'^^f{x) (17) 

We also set 

so that at{x) = e"*'^'/^. Let us now give explicitly the terms in (llSp . 

The source terms (that is, the terms in the left hand side of the Bethe equations) in the case 
of alternating 3 and 3 ([6]) for n wide strings of type I or type II are of the form a2n-20t2n which 
gives, after taking the log and differentiating, source terms of the form 



•Sj7,n = «2n + a2n-2, n>l (19) 



with the convention that oq = 0, and that wide strings with n = 1 are just real strings. Similarly, 
the source terms for n strange strings of type I or type II are still complex numbers of modulus 
one, e2n-ie2n+i) and give after taking the logarithm and differentiating, real source terms of the 
form 

^l,n = Sh,n = a2n-l + a2n+l, n>l (20) 

The kernel of interaction between n wide strings of type I and p wide strings of type I reads 

«2p-2na2p-2n+4 • • • 0^2p+2ra-4 (21) 

while the interaction between n wide strings of type I and p wide strings of type II is instead 

a2p+2n-2a2p-2n-6 • • • «2p-2n+2 (22) 

and there is a symmetry between type I and type II strings. 
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Taking the logarithm and differentiating gives the Fourier transforms of the functions $ in 
(seeds]) 

iw;w _ iw;w _ Smh nx rp_n)\x\ 

^I,n;I,p - ^I,p;I,n " ^ ginh X - 

^>/,;T/7,P = ^7i:;;/,n = -^'-^^-^'^^ - e-(^+"-^)l^l P > n (23) 

^Ik!/,, + Kn.n,p = -^2n-l,2p-l + 5n,p (24) 



Note that 
where 



^r,. ^2^5^J^sinh[rx/2]e-^l^l/2, s > r, = 
smhx/2 



(25) 



and —Ars + Srs is the Fourier transform of the standard kernel describing interaction between 
r and s strings in the standard Bethe ansatz for XXX chain. Of course, this result should not 
be a surprise: concatenating u and 7 variables transforms a n wide string of type I or II into a 
2n — 1 string for the XXX equations. Meanwhile, the source terms also match, for the chain acts 
on spin one representations. 

We now turn to strange strings. Some of the kernels then will not be real, but after concate- 
nation of the solutions, they will correspond to strings with even labels in the XXX model, a 
feature already seen for the source terms. For instance the interaction of type I n strange strings 
with type I p wide strings corresponds to the kernel 

[u + i{p-n- I)] [u + i{p-n+ 1)]^ ...[u + i{p + n- + i{p + n - |)] ^^^^ 

[u — i{p + n — |)]2 . . . [u — i{p — n + 

Meanwhile the interaction of type I n strange strings with type II p wide strings is exactly 
the inverse of the complex conjugate of this quantity. After taking the usual logarithm and 
differentiating one finds the Fourier transforms of the functions <I> (jlSp 

^/,n;/,p = \ [-Mn,2p-l - Sign {x)B2n,2p-l] 

^'l%;n,p = \ [-Mn,2p-l + sign {x)B2n,2p-l] (27) 



where 



2sinhx/2 

Observe the sum rule 



i?r-,s = ^^j^sinh(rx/2)e-^l"l/2, s > r > 1, 5.,. = 5.,. (28) 



^li,p + ^iln,r> = -^2n,2p-i (29) 
The interaction with real strings is a bit different, as 

i?i,, = 2sinhMe-N/2 (30) 
Finally, we need the interaction of strange strings among themshelves. One finds 

is;s — _C -I- A 
^I,n;I,p ~ ^i^^P ^ "n,p 

^TnJI,P = -Dn,p (31) 
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with 



2 cosh X 



sinhx 

sinh X 
Note that, once again, 



Cp,n = —r-, sinh nxe pI""! , p>n, Cp,„ = C„ „ 

smhx 

Dp^n = . \ sinhnxe"^!""!, p > n, Z)p,„ = Dn^p (32) 



We now write symboUcally the system of equations as in (I15p but allowing explicitely for the 
two types of solutions: 

Phi + P\,i = + * Pl,j + * Pll,j 

j 

Pll,i + Pn,i = sn,i + * Pl,j + * Pll,j (34) 



and the labels i, j are short hands for n,p and s,w. The energy is given by (we often do not write 
the integration variable) 

e/,i pi^i + en,i Pii,i (35) 



E 

T 



where ^ = efj^ = —2sf^ = —2sfj^, and similarly for the strange strings. Note that although 
the Bethe solutions for strange strings are not invariant under complex conjugation, the symmetry 
between the two types of roots in the expression of the energy guarantees that it is real. Note 
we have adjusted the sound velocity (using the results of [8j for low energy excitations) so that 
the theory is isotropic in the continuum limit. 

We must now pause to discuss the effects of the strange strings more carefully. Since some 
of the corresponding bare scatterings are not actually pure phases, the associated kernels are not 
real (this translates into some of the Fourier transforms having an even and an odd part under 
change of sign of the variable x; see (j27p for instance). On the other hand, the densities are 
defined as variations of integers with respect to a real variable (the center of the strings), and 
thus are real. The equations ()34p can only make sense if the right hand sides are real as well. It 
seems reasonable to solve this constraint by demanding that the densities of type / and type // 
n strange strings are equal (as functions of the rapidities). 

The densities being real, the entropy has the usual expression in the TBA: 



S 
L 



j YipLi + Pii) HPl,i + Pii) - Pl,i In Pl,i - Pi^ In P^j +1^11 (36) 



Note here that one has to be careful to not constrain the individual type / and type II strings 
to have the same centers, as this would reduce the entropy by a factor two, and give rise to the 
wrong result in the end. Rather, we use the expressions (j34|[7l [36j) and minimize the free energy 
F = U — TS as a function of the densities pi^i,pn,i subject to the constraint that the densities 
of strange strings be equal. This is in fact a transparent subtlety, if the naive extremum of F 
without any constraint turns out to be symmetric under / •(->• //. 



11 



Introducing pseudoenergies through p\l Pn = e^"/^, we find the (formal) TBA equations 
= e/,i - T ^ ^i^i.jj * In (l + e-^^.^/^) - T ^ ^>,,,,,,,- * In (l + e-^">^/^) 

eii,i = en,i -TJ2 ^i,r,ii,j * In (l + e"^^-^/^) - T J]] * In (l + e"^"-^/^) (37) 

while 

T = -^Y. f [^^'* 1^ + e"'^-'/^) + In (l + e-^"-/^) ] (38) 
Demanding e/^j = en,i = at equilibrium, we get a new system which reads 

ei = ei-Y^ T{^i,r,i,j + * In (l + e"^^/^) (39) 



together with 

F 
T 



-2TY, j Si\ii(l + e-'^/^^ (40) 



Using the sum rules mentioned previously, it is clear that the resulting TBA equations coincide 
with those for the spin one XXX chain, while the free energy is twice the free energy of that 
chain, leading to a central charge c = 2x|=3in the UV. 

Let us be more explicit on these points. By simple manipulations ones finds the universal 
form of dSD 

2 cosh I {pj^i + pii^i) = 25i,2 + + P^li+iP^ii,i-i + (41) 

and similarly the universal form of ([3 



4s 

e/,i + e//,j = 5j,2 + Ts-k 

TT 



In i^l + e^'-'-^/^j + In (^1 + e'-^'^+i/^ 
+ ln fl + e^".»-i/^^ + In fl + e'"-'+^/'^ 



(42) 



with s' = 2cosha72 ' ''^^ Convenient to represent these equations by a TBA diagram as shown on 
the figure [9l Dots stand for wide strings and squares for strange strings. Note how they nicely 
get organized into a pattern exactly identical to the one of the XXX chain. At T = the ground 
state is obtained by filling up states corresponding to the black squares. Of course the presence of 
two lines of nodes corresponds to the algebra being of rank two; moreover the symmetry between 
the two lines arises from our symmetric choice of roots. 

It is a simple matter to extend these results to a chain with alternating fully symmetric 
and antisymmetric representations, represented by Young diagrams with p boxes . The TBA in 
universal form is represented by diagrams identical to the foregoing ones, only the massive nodes 
are the 2p*'' ones. The central charge in the UV is c = 

It is also easy to extend the results to the case of the fundamental representation. The TBA 
now has source terms on the first nodes, as represented on figure [TUl The UV central charge is 
c = 2. 

The TBA for the spin chain confirms identifications of central charges carried out by other 
means in [8], [18]. Moreover, it constitutes a crucial check of completeness of the solutions of 
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Figure 9: TBA diagram for the s/(2/l) chain with alternating 3 and 3 representations. 




Figure 10: TBA diagram for the OSP{2/2) chain (four dimensional fundamental representation). 

the Bethe ansatz equations, and we can now use our strings classification to tackle the more 

interesting problem of scattering in osp{2/2) intcgrable quantum field theories. 

Nevertheless, it is fair to remind the reader that the TBA is rather insensitive to details of 
the diagonalization. Consider for instance the well known SU{2) spin 1/2 XXX chain. The total 

number of states with 5*^ = is ^s^=o = ^ ^ while the number of highest weight states with 
_7 = is flj=o = ( ^l'^ ] ~ = ( ^'^^r ^ ] • However, if we consider the entropies associated 



^ L J \ L ^ 
with these two numbers we have 

L + 2 

In r2j=o = In ^s''=o ~ In — ^ — (43) 

so the difference is of order and goes to zero as L — )■ oo. In other words, thermodynamics 
based on a scheme that would not provide all eigenstates but only the highest weight states say 
would still reproduce the expected results. In particular, even if our string analysis turns out to 
be wrong or incomplete for states in indecomposable representations, it is still very conceivable 
- as we indeed observe - that this effect is negligible in the thermodynamic limit. 

A related comment concerns the reality condition we imposed on the densities of the two types 
of strings. Since evidence is strong that the minimum of the free energy functional is symmetric 
under the exchange of the two types of strings, our condition led us to what should be the correct 
result. It is not clear however how to describe solutions where there is a dissymmetry between 
type I and type II strange strings. Presumably, the string hypothesis breaks down in this case. 
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4 TBA for OSp(2|2) scattering: the fundamental representation 
case and the GN model 

We now consider the different but related problem of writing up the TBA for the OSp(2\2) GN 
model. The action is 



(fx 



■ 2 



(44) 

where the if^ are Majorana fermions of conformal weight 1/2, and the (3^ are bosonic ghosts of 
weight 1/2 as well. The theory at g = can be identified with the OSp{2\2) WZW model at 
level k = 1 (part of the literature denotes this level k = —1/2 instead, focussing on the sub 
SU(2)). The lowest order beta function is j3{g) oc —g"^ and for 5 < the theory is massive B The 
corresponding S matrix has been proposed by Bassi and Leclair in [12]. It involves, as is usual 
for 0{N) GN models at large but not, a priori, at small N , particles of mass M in the four 
dimensional fundamental representation, with scattering matrix of the form S{9) = Z[6)R^^^^{9) 
where Z is a normalization factor, and -Rfm^jj is the solution of the Yang Baxter equation for 
the fundamental representation, normalized such that the scattering kernel for particles in the 
highest weight state b = 0,S^ = 1/2 is equal to Z (that is the matrix element in R is equal to 
unity) . The R matrix is the same as the matrix defining the Boltzmann weights of the integrable 
spin chain, in the fundamental representation. The minimal Z factor (in general we denote by 
Z factor the element of the S matrix corresponding to diagonal scattering of two particles of the 
same type - no ambiguity will arise in what follows) given in [12j using standard crossing and 
unitarity arguments turns out to be the square of the similar factor for the SU (2)i WZW theory, 
or sine-Gordon model at = Svr, 



r(i + f)r(^-if) 
r(i-|f)rQ + i^) 



(45) 



Doubts were raised in [Ij about this S matrix, on the grounds that, for OSp{2\2) as for the 
usual 0{N) GN models at small N, the fundamental particles might well be kinks, while the 
fundamental fermions (and bosons in the super case) are unstable. It was suggested that maybe 
the S matrix of [12j in fact describes the sigma model OSp{2\2)/OSp{l/2) instead. We will, in 
what follows, find evidence that this expectation is not correct, and that the S matrix of |12] 
does indeed describe the GN model. Of course, writing the TBA for the scattering theory is 
not the same as writing it for the integrable spin chain, as the number of particles is allowed 
to fluctuate. The two problems are however closely related, as is well known (see |19] for a 
pedagogical discussion) 

The problem of quantizing the particles turns, in the continuous limit, into the basic equation 

2Tr(a + a^) = M cosh9 + --^InZ a + -^InA (46) 

I dO I dO 

Here A is the eigenvalue of the monodromy matrix, ie of the matrix describing the process of 
passing one particle through a system with a certain number of particles already. This eigenvalue 
(properly normalized, so that it is equal to unity if all particles are of the same type, corresponding 



*We use Majorana instead of complex fermions. Some of the signs are therefore switched when compared to 
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to 'all spins up' in the auxiliary problem) is found itself by another, 'auxiliary' Bethe ansatz, which 
coincides with the one used to diagonalize the spin chain. We need the general formula for the 
eigenvalue picked up when passing a particle of rapidity v through a set of particles with rapidities 

{va} 

(one checks that its logarithmic derivative gives the energy for the hamiltonian in the previous 
section), where the n,7 are related to the {va} through the Bethe equations 



M 

A J ^ ' 13= 



N 

jj 7.-^. + V2 jjT,,-^^ « = 1,...,M (48) 

Note that the same set of rapidities {va} is used in both equations for the u as well as the 7's. 
Here we have used the lattice rapidities for the auxiliary problem to match with the previous 
section. They are related to the relativistic ones by ^ = ttu. We define the Fourier transformation 
for functions of 9 as 

fix) = [ dOe^'-l^fie) (49) 



The Bethe equations for the auxiliary system (|48p admit the same clasification of roots as the 
one for the lattice model, and lead therefore to the following system 

Pl,n + Pl,n = a2n- 1*^ + 2^ */,n;/,p * Pl,V + ^I,n;II,p * Pll,p 

P 

+ ^Z-,i,p^pip + ^Z-j,P^ph,p 

Pin + Pin = <^2n^ a + Y,^7Zi,P*Pf,p + ^Zn,P*Pn,P 

p 

+ '^l%;I,p*Pf,p + '^Z;I,p-*^Pll,p (50) 



and a similar equation for type II roots (the kernels $ have been defined in the previous sections) . 
As for the a equations they are H 

a + a'' = ^ cosh e + ^^lnZ-ka-Y, a2k-i * {p7,k + Pn,k) 



27r 2i7r dO 

k>l 



^a2k*{pj,k + Pii,k) (51) 



k>l 

where we now define 

«t(^) = -:|ln|^^ (52) 

^ ^ 2Trde e-iTrt/2 ^ ' 

and all the densities are now in terms of the (physical or auxiliary) rapidities. 



^The mass parameter M should not be confused with the number of roots of type 7. 
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We write symbolically our system as 

Pi,i + Pti = ^i^i;^ * + XI ^I'^-'^'j * Pi'i + ^i,i;ii,3 * Pn,j 

j 

Pll,i + Pn,i = * 0- + X] ^Il'i'Ij * Plj + * PIIJ 

3 

Mcosh.6 ^ \ - , ^ , 

+ fj'' = + ^a,a * 0" + 2^ * + * Pll,i (53) 

i 

and the labels i, j are short hands for n, p and s, w. We analyze then the thermodynamics of this 
system, with the energy 

E 



L 



j (7 {6) M cosh Ode (54) 
Introducing pseudoenergies through / p = e*^"/"^, /a = e^*^/"^, we find the TBA equations 
e = M cosh 9 - T^^,^ * In (^1 + e"'-/^) " ^ * 1^ (^1 + e"'^-*/^^ 

- T$7/,i;^*ln(l + e-^"-/^^ 
e/,i = -r$,;,,, * In (l + e-^-/^) - T * In (l + e-'''^'^" 

3 

- In (l + e-^".^/^) 

3 

en,^ = -T^a;ii,i * In (l + e"^-/^) - T «>j,,;77,i * In (l + e"^^-^/^) 



3 



while 



T J2 ^n,r,ii,i * In (l + e-^"-^ /^) (55) 



^ = -r/^^^|^ln(l + e--/^)d0 (56) 



Observing now that = = ^(T,i, we expect that e/^j = e//^j = at equilibrium, and 

thus get a new system 

= - 5^ T{^i,r,i,i + ^ii,3 j,i) * In (l + e"^^/'^) - T^^.i * In (l + e"^-^/^) 
e^ = M cosh - r$^,^ ★ In (^1 + e"'"/^) - 2^ ^ ★ In (^1 + e"''/^) (57) 

i 

Remarkably, this TBA is identical with the TBA one would write for the SU{2) PCM model 
|20j . There, the scattering matrix is the product of two SU{2) isotropic sine-Gordon scattering 
matrices, and the strings involved in diagonalizing each of the two scattering matrices behave 
like our type / and type II solutions. 

Identifying the nodes of type I and type II and using the fact that the pseudoenergies are 
equal in equilibrium gives the equivalent diagrams shown in 1121 
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Figure 11: TBA diagram for the OSP{2/2) theory with particles in the fundamental representa- 
tion. 



■o 



o — o — • — o — o 



Figure 12: Equivalent TBA diagrams for the OSP{2/2) theory with particles in the fundamental 
representation, showing mapping to the SU (2) PCM model. 
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We conclude that cuv = 3, for our scattering theory. This is in agreement with the effective 
central charge for the GN model in the sector with anti-periodic fermions (NS), and provides 
some justification for the S matrix of [12j . 

If we now accept that the S matrix in [12] is the correct one, it follows from our calculation 
that the free energy of the OSP{2/2) GN model in the NS sector is exactly equal to the one of 
the SU{2) PCM model. Now note the latter is the same as the SO{A) / SO{?>) sigma model: we 
thus have established part of the equivalences expected from the introduction. 

Our observation deserves several remarks. The first is, that the torus partition function of 
the free GN model with periodic boundary conditions for the fermions is exactly equal to unity. 
This results from the cancellation between the fermionic and bosonic contributions (each equal to 
Det A^^, A the Laplacian), and is in fact what motivates the appearance of the OSp{2\2) symme- 
try in the SUSY approach to disordered systems. One can easily show that this supersymmetry 
is not broken by the introduction of a current-current perturbation. The partition function with 
periodic BC should thus remain unity, which implies in particular that the hulk term in the free 
energy - which does not depend on the boundary conditions - should vanish exactly. Turning now 
to the non trivial sector with antiperiodic boundary conditions for the fermions, the free energy 
now becomes the non trivial object we have studied in this section, but its bulk part should still 
vanish exactly. Since we have argued the free energy for the OSp{2\2) GN model should coincide 
with the one of the SU{2) PCM, this implies the bulk term must vanish there. To calculate this 
bulk term, we use the fact that SU(2) PCM coincides with the 0(4) sigma model. The bulk term 
for 0{N) sigma models follows from the general results in [21] , [12] ; [23] : 

— — = cot — 58 

L 8 N-2 ^ ' 

where M is the lowest mass. We conclude that for the SU (2) PCM, for which = 4: 



L 

in agreement with our expectation for OSp{2\2). 



X (59) 



5 The OSp(2|2) supersphere sigma model TBA. 

If the S matrix proposed in |12j is for the GN model indeed, the question is now open as to 
what the S matirx of the sigma model might be. By analogy with what happens for 0{N), 
N > 0, and after switching the roles of GN and cjM, one would expect that the S matrix for the 
OSp{2\2) sigma model should involve spinor representations. This could be the three dimensional 
representations (fundamental of 5'L(2|1) and its conjugate), or maybe some infinite dimensional 
representations [2l]. Our naive attempts in this direction have not been successful however. 

An alternative approach consists in starting from lattice regularizations. It was indeed argued 
in pp that the continuum limit of the integrable staggered OSp{2\2) chain built on the fundamental 
was described by the OSp{2\2)/OSp{l\2) sigma model in the massive regime. The staggering is 
similar to what has been studied in particular under the name of lattice light cone regularization 
for a large class of perturbed CFTs and sigma models [25j [26j . For a recent reference in a very 
similar case, see [27] . 

Recall that, for the integrable chain based on the four dimensional representation, the bare 
source term involves a factor | instead of i, and as a result, the source term for the wide and 
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strange strings are different: 



S/,n — ^11,71 — 0.2n-l 

S^n = Sji^n = a2n (60) 

(so in general we can write that the source term of a string with p roots is Op). 

What we have to do, following the general ideas of [lj,[27j is to consider instead of the 
homogeneous chain, the chain arising from staggering spectral parameters. The relevant equations 
now read 

Uj + A/2 + i/2\^^^ fuj - A/2 + i/2\^/^ _ tt ""i " 7/? + ^ 



Uj + A/2 -i/2 J \Uj - A/2 -i/2 J ""i - 7/3 



+ A/2 + V2y/^ - A/2 + ^/2y/^ ^f.,.-u, + ^ 
\-fa + A/2-i/2j V7a-A/2-V2y 7a " «fc - i ' 



k=l 

with the quantum numbers B = ^"^^'^ and = L — The energy reads then 

/f 1 1 f 1 

(^fc - A/2)2 + 1/4 + (ufc + A/2)2 + 1/4 + (7^3 - A/2)2 + 1/4 + (7^3 + A/2)2 + 1/4 

(62) 

The ground state of the system was discussed in [HI [18] and was shown to consist of filled up sea 
of real roots (ie n = 1 wide strings) of type I and type II. Excitations can easily be shown (see 
more details below) to be massive, with a mass scale Ma oc e"'^^/^, with a the lattice spacing, 
so the continuum limit is obtained as usual with a — ?■ 0, A — )• c«. From the general discussion 
in the second section, it follows that the ground state energy is given by twice the ground state 
energy of the staggered XXX chain. But for the latter, the ground state energy is well known to 
reproduce, in the continuum limit where A — t- 00, the ground state energy of the SU{2)i WZW 
model perturbed by the current current interaction, or, equivalently, of the sine-Gordon model 
at = Sir. 

It is also a well known fact that the 5*0(4) GN model is equivalent to two decoupled sine 
Gordon models at their marginal point [28] . We thus have obtained our other claim, namely that 
the free energy of the OSp{2\2) supersphere sigma model with antiperiodic boundary conditions 
coincides with the free energy of the 50(4) GN model. 

We thus have coincidence of free energies as 

OSP{2/2) GN o 50(4)/50(3) sigma Model = SU{2) x SU{2) PCM 
OSP{2/2)/OSP{l/2) supersphere sigma Model o 50(4) GN = SU{2) x SU{2) GN (63) 

We note however that the relationship does not seem to extend to the full spectra of the 
theories. We haven't been able to find an SU{2) x SU{2) symmetry within either of the OSP{2/2) 
models for instance. The UV limits of the flows are quite different as well. For instance the UV 
limit of the 50(4) GN is compact, while the UV limit of the OSP{2/2) supersphere sigma model 
has a non compact direction. 

Now we can get back to the question of what the S matrix might be. For ordinary models, this 
can be rather easily determined using the Bethe equations from the lattice regularization. The 
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idea is to rewrite the continuous Bethe ansatz equations in terms of excitations over the vacuum 
p9] . ie by putting on the right hand side, not densities of n = 1 wide strings, but densities of 



n 



1 holes in the wide strings ground state distribution, that is replace pi^i and pn,i by Pi \ 
andp^i|\. 

After laborious but straightforward calculations, one ends up with the following equations for 
the densities of real particles (we denote aj = pf-^ and similarly for type II) 



coshAx/2 e l^'l ^ 1 ^fe (fc_i)|^|^ 

+ = 2cosh(x/2) - 2^hIhN ' + 2^kdrN^^^ " P''' 

^ {k-i/2)\x\ ( 1 + sign (x) 1 - sign (x) \ 

k>\ ^ ^ 

, ^fe cosh Ax/2 e-l^l ^ 1 ^fe 

^ ' ^ II II k>2 

- E f^^^T^a. + i^^T^ftO (64) 



k>l 
COS 

2cosh(x/2) 2sinh|x| '''' ' 2sinh|x| 

^ ' ^ II II k>2 

' 1 + sign (x) ^ 1 - sign (x) , 

^ V - " ' 

fc>l ^ ' 

For the other solutions, one observes that the bare source terms disappear, and are replaced by 
source terms which are convolutions involving the densities and <Tjj. As for the interactions 
between different types of roots, one finds that a shift of one unit in the length is taking place. 
That is, what was before an n (n > 1) wide root of type I or II now behaves as an n — 1 strange 
root of type I or II, while what was before a strange n root of type I or II now behaves as an n 
wide root of type I or II. In other words, we have the following system 



Pl,k+Pl,k = +l^^^I,k;I,l-*' Pl,l + ■■■ 

l>2 

Pii,k + Pii,k = e ^ ^ajj + 2^D<^jjf^.jii<pj^i + ... 

l>2 

Pi,k + Pi,k = e ' / ^1 I I aj + ajj I + + 2^ ^*/,fc;/,z * Pi,i + ■■■ 

^ ^ l>2 

.-^fA ,-(fc~i/2)k| 1 + sign(x) l-sign(x) fe^^ , ^ T^<f.s,u, 
Pii,k + Pii,k = + + + Z^D<^iLkj,i-*'Pi,i (65) 

^ ^ l>2 

where each time the dots indicate the obvious sum over the remaining terms (/, s), {II, w), {II, s), 
and the dressed kernels are obtained by shifts 

DK:Z,P = ^rn-i;/,P-i • • • (66) 

The new system bears a lot of resemblance with the old one. In fact, up to a relabelling of the 
roots, it is identical for what concerns the interactions between the densities p of pseudoparticles. 

This is fully expected, since the dynamics of the p excitations is determined entirely by the 
symmetries of the model - here, osp{2\2) super Lie algebra. 

Now for the interaction between the sigma densities. In the case of ordinary Lie algebras - 
say sl{3), of rank two like our present TBA - the holes in aiji distributions would be associated 
with particles living in the representations corresponding to the first top or bottom node of 
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the TBA diagram. Oddly enough, while this would be fine if we were considering the second 
nodes on the diagram (which correspond, as we have seen earlier, to 3,3), there is no osp{2\2) 
representation corresponding to the first nodes individually! A quick calculation shows that these 
representations should have dimension s/i = 2 indeed, and this is not possible in osp{2\2). This 
might indicate that the full osp{2\2) symmetry is not present in the S matrix, but maybe only a 
sub sl{2) X .sZ(2). 

However, it should also be clear that the kernels for I — I and I — II interactions are also 
quite ill behaved, as their Fourier transforms both diverge at the origin: 



I d ^ ^ 1 d , ^ If, e-l^l 

I dd t dd TT J 2 smh | 

1 d , ^ lci.„ If, 1 



. In Zjji = ^ ;^ In Zjj,! = - / dx e-'^^l^ (67) 

I dQ I dQ IT J 2smh|a:| 

(here, 9 is the usual physical rapidity, obtained from the lattice rapidity through a rescaling). 
We believe this divergence is related with the non-compactedness of the target space in the UV 
limit; this will be discussed in more details in the next section. 

Meanwhile, it is intriguing to observe that, if we were to define new excitations by binding a 
hole in the distribution of real wide roots of type I and a hole in the distribution of wide roots 
of type II at the same rapidity, the Fourier transform of the derivative of the logarithm of the 
dressed scattering kernel of these new excitations with themshelves would be 

/ 1 e-l^l \ 

^nii,i+ii = 2 (^^^7^ - ^^T^ J = ^^^^ (68) 

Similarly the kernel of such an excitation with a wide string of type I or type II is given by 
e^C^-i)!^! while the kernel of such an excitation with a strange string of type I or II is given 
by e('=-V2)l^l. Hence, if we forced holes in the I and II distributions to be bound, the following 
system of equations deduced from the lattice model would now coincide with the TBA of the 
previous section provided we identify 

27r$,,, = , In Z = - / -^^^e-'^-l-dx (69) 
I dd TT cosh(a;/2) 

Integrating this quantity we find 



r(i-if)rG + |f) 



(70) 



which in fact is the result for the GN model discussed above. In other words, pairs of holes in 
the I and II distributions at identical rapidities in the sigma model scatter exactly with the GN 
S matrix! There is however no obvious reason why these holes should be paired (and the sigma 
model does not coincide with the GN model anyhow), so the sigma model S matrix appears as 
some sort of 'split' GN S matrix. More work on the bound states and crossing unitarity would 
be necessary to give sense to this idea. 

We note here that in the case N > 0, the generic sigma model S matrix gives the correct 
scattering theory for all values of N, while the generic GN S matrix works only for large enough 
N; at small enough A'', the fundamental fermions are unstable, and the correct S matrix involves 
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instead kinks in the spinor representations. It seems that for N < 0, the roles of the GN and 
sigma model are switched again: the generic S matrix for GN seems to still work for = 0, 
while for the sigma model at = 0, some sort of 'spinor' S matrix seems necessary. The spinor 
representations of OSp{2\2) are however infinite dimensional, and it is not clear whether they 
would allow for a meaningful description of the scattering in the present case. 

We note finally that the ground state energy of the OSp{2\2) sigma model can be calculated 
using the general result (I58p with = and one finds again 

This agrees with the ground state energy for the WZW model perturbed by current current 
interaction [30] . 



6 Do the massive flows exist? 

A important remark is that, according to the discussion in |31] . the OSp{2\2) WZW model 
perturbed by a current current interaction in the massive direction in fact does not exist as a 
field theory. More precisely, it is established in this reference that the functional integral is not 
well defined, the topological term leading to strong divergences. 

The same holds true, in fact, for the sigma model, as was briefly discussed in [32\ . The action 
for the supersphere OSP{2/2) /OSP(l/2) sigma model is obtained as follows. Parameterize the 
5^'^ supersphere as 

= COS0 (1 — rjir]2) 
4)2 = sin0(l - ?7ir?2) (72) 

such that (/>f + 02 + 2r/i?72 = 1- The action is then 

5* = ^ y d^x {df,(t}f (1 - 2r?i?72) + 2df,r]id^r]2 - ^mmdiiVidi^m (73) 

where (j) is compactified, = </> + 2tt, and conventions such that the Boltzmann weight is . 
The spontaneously broken symmetry phase, with the theory being free in the UV and massless 
in the IR, corresponds to a positive coupling constant g, since for the 0{N) model, the beta 
function is 

m = {N-2)g^ (74) 



dl 

and the same holds for 0Sp{n\2'm) sigma models with N = n—2m. When g is positive, the theory 
flows to weak coupling, and critical properties are described by the 5 — )• limit. A rescaling and 
relabeling brings the action into the form 

S = j d^x {d^4>f (1 - 2gr]ir]2) + 2d^rjidf,r]2 - ^griimdfxriidiMm (75) 

with now (f) = As(7— )-0, we thus obtain a symplectic fermion and a non-compact boson. 

Note that the presence of a continuous spectrum in this osp{2/2) version of the problem has 
a physical origin in the fact that the symmetry is spontaneously broken, and thus correlation 
functions of order parameters have no algebraic decay (though they have non trivial logarithmic 
behaviour) . 
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Now the massive region meanwhile corresponds to the theory flowing to strong couphng, and 
g negative. Except in the OSp{l\2) case, where there is a single bosonic field, and the theory can 
be reformulated entirely in terms of the symplectic fermions, in the 05^(212) case in particular, 
a negative coupling constant leads to a term exp [+ f cPx{d^(p)'^~\ in the Boltzmann weight, and 
thus to a theory that is not well defined from the path integral point of view. We note that 
trying to define theories in a region where the path integral is not defined is a task that has 
been considered in string theory for some time. An example is the so called "time like Liouville 
theory" , discussed for instance in [33] and [31] . 

So what happens at g negative and the "massive" region is not entirely clear. The model 
defined by the action ([73j) is unstable, a fact that probably manifests itself by a spectrum that 
becomes unbounded from below, and hence level crossing right at g = 0. 

The point can be illustrated more concretely in the case of the 05*^(412) GN model. In the 
papers |5] and [6] it was found that, for the sign of the 0{N) GN coupling which corresponds to 
the relevant direction in the case N > 2, the model is gapless and exhibits a line of fixed points 
with central charge c = 1. Meanwhile, for the opposite sign of the coupling, something very 
different occurs. This can best be seen by contemplating the formula giving the exact values of 
the conformal weights and thus the scaled gaps on the critical line (for the theory with Neumann 
boundary conditions), as a function of the Young diagrams for the 05^(412) representations. 
Writing the weights at the free point as 

n? Tin b — 2 
^free = Y + f + 

2 

^free = Y' ^ = 0' ^ C^S) 
for a Young diagram with shape A = n\n2l^^^ , the exact formula [5l|6] is 

h = h^j-QQ + gC (77) 
where C is the Casimir of the representation 

^^(m+^W)!^(m^^+i)^_^^_^^2 (78) 

While for g > the ground state remains in the identity sector, for g < 0, the weights of the 
fully antisymmetric representations 1^ (for instance) become arbitrarily large and negative for 
p = 0{1/ yfg) as soon as g is turned on: an infinity of crossings occurs, and the new ground state 
has nothing to do with the one at 5 < 0. 

On the other hand, we could also restrict to the 0(2) subsector of the theory. There, the 
partition function is the one of the free boson with Neumann boundary conditions, and nothing 
happens - since all the representations which could cross the ground state are eliminated when 
going from OSp{A\2) to the sub 0{2). The properties of the ground state, analytically continued 
in this phase, are well behaved, but describe a subset of the theory which is very high in energy. 
The equivalent observation in the OSp{2\2) GN model would correspond to the sector with 
periodic boundary conditions, where the partition function remains equal to unity and c = 0; 
this sector might be very high in energy with respect to the true ground state in the would be 
massive phase. 

The 05^(212) sigma model can be studied somewhat in the same spirit by using perturba- 
tion theory. Near the weak coupling fixed point, the massive direction formally corresponds to 
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negative g, and, if we for a minute do not worry about the meaning of the path integral, we 
can calculate the evolution of the hamiltonian eigenvalues perturbatively. Symmetric representa- 
tions are particularly easy to handle, since they can also be obtained through the minisuperspace 
analysis. The eigenvalues are simply proportional to those of the Laplacian on the supersphere: 



and of the form E = gn{n — 2), n an integer. They all collapse at 5 = 0, and there is an infinity 
of level crossings at that point when going from g positive to g negative. 

We note that, if this scenario is correct, the point g = must be considered as a 'first order 
critical point". This is a point where the theory is conformal invariant, but which is at the same 
time a point where level crossing occurs, such that the analytic continuation of the ground state 
becomes a very highly excited state in the low temperature phase. In general, first order phase 
transitions are not associated with conformal field theories, but several examples of first order 
critical points are known. These involve crossings in finite size, and are only possible because of 
the non unitarity. Maybe the best studied example corresponds to the thermal perturbation of 
the antiferromagnetic Potts model on the square lattice, as discussed in [18j. 

Meanwhile, we have a lattice regularization (the staggered spin chain) and a TBA that seem 
perfectly stable. In particular, we have not found, by exploring solutions which are dissymmetric 
in the two types of roots I and II, any indication that the ground state of the chain is not given 
by the foregoing results. It looks thus as if our results describe some sort of continuation of the 
unstable theory, which now behaves as an (almost) ordinary massive theory. It would of course 
be interesting to explore this issue in more details. 

7 The q-deformed case 

We have so far discussed only the case of twisted boundary conditions for the fermions, corre- 
sponding to the "effective" central charge of our non unitary theories. To get the true central 
charge, one must consider a twisted TBA, by inserting (complex) chemical potentials for the 
kinks in the scattering. For the OSp{2\2) GN model, we do not expect that much interesting will 
remain. The central charge will be zero in the UV, it is zero in the IR for the massive flow, and 
probably remains zero all the way. 

The situation is more interesting for the supersphere sigma model, since now c = —1 in the 
UV. To get the calculation of the true central charge under control, it is convenient to consider 
the q deformed version of our problem. 

In terms of lattice model, the q-deformation is obtained by taking the four dimensional solu- 
tion of the YB equation based on UqOSp{2\2) 

Setting q = e*^ we have the Bethe equations for a chain with L sites 



(the isotropic limit is recovered by setting Uj = 2^Uj,r^ = 2^7^ and ^ — > 0. The energy of the 



H = gAs 



'-1 



(79) 





(80) 
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chain is meanwhile 

E I ^^^^ ST^ si^C I /g-|^\ 

1 ^ cosh Uj — cos ^ ^ cosh — cos ^ I 

In the antiferromagnetic regime ^ G [0, |], the ground state is obtained by fiUing up the seas yi, yu 
with real roots. With anti-periodic boundary conditions for the fermionic degrees of freedom (the 
NS sector), the central charge reads Cgg = 2. 

The most natural is to observe that the XXZ sub component of the spectrum (that is, the 
component given by yi = yu) can be twisted to obtain a c < 1 theory described by a Coulomb 
gas with parameter ao- To do so we impose boundary conditions twisted by e^*'^'^ on the spins, 
^ = f resulting in the central charge for the XXZ chain 

^(«/^'' (82) 



1 - 

We now parametrize 

1 

(83) 



(usually one would set ^/tt = l/(m + 1)) so the central charge of the OSp chain becomes 

6 



^twistl - 2 



1 - 



(2I + 2)(2I + 3)_ 



(84) 



Meanwhile, we recall that the effective central charges (ie, describing the ground state scaling in 
the NS sector) of the 05^(212) and OSp{l\2) WZW theories at level I read 

OSp{2\2)i; c = 0; c^S = 

91 8/ 

OSpmr, -aT3^ 'eff=^ (85) 

It is easy to check that 

6/ 8/ 

ctwisti - [qri - ^ (86) 

This suggests that the continuum limit of the spin chain is described, for I integer, by coset 
models OSp{2\2)i/OSp{l\2)i and suitable interpolations in between. 

When one goes from the XXZ to the minimal models, a set of electric charges have to be 
introduced, which are multiple of the fundamental charge eo = ^. The particular choice e = 
(Z + l)eo gives the conformal weight, with respect to the central charge of the minimal model, 

/..= 



4(2Z + 2)(2Z + 3) 
This leads to an effective central charge for the OSp theory 

21 

ctwist2 = ctwistl - 24 X 2 X /i = - (87) 

which is the true central charge of the coset model. The phase shift to obtain this central charge 
is thus cxp[2i7r(Z + 1)/ {21 + 3)] ; it becomes a pure e*'^ in the I oo limit, where the central charge 
c^^jg^2 becomes equal to —1. 
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A way to obtain the flow of the central charge in the supersphere sigma model is thus to take 
the flow of the running dimension in the minimal models M2i+2,2l+3 iii the limit ^ — )• oo. 

Now this is particularly simple. Indeed, we observe that ht = and the running 

dimension of operators in minimal models is the easiest to obtain, since it follows from a 
simple twist on the basic sine-Gordon equations. This can be done in the context of the TBA, 
as discussed in |36]. In the notations of this paper, one must set 7 = |, so the soliton fugacities 
are exp[ibi7r(2/ + 3)/2]. One can also use the DDV formalism, and for instance the formulas in 
(37] . The basic equations reads 

f{0) = imR sinh 9 + ia + j -^{O - e')ln (l + e-^^^'^) dO' + j <^>{e - 9') In (l + e'-^^^'^) d9' (88) 

where ^> is the SG kernel for the corresponding value of the sine-Gordon coupling, = 1 — 51+3 ~ 
11^, Ci and C2 are contours slightly below and above the real axis, and the effective central 
charge follows from 



^ 2>imR 
c{mR, a) = — ^ — 



s\n\ie\n{l + ef^^^^d9 + j s{nh.9\n{l + 



Ci JC2 

I 



d0 (89) 



In the present case we need a = f exactly. The flow now should go from c = — 21^:3 in the UV to 
c = in the IR. In the limit / — t- 00, this becomes a flow from c = — 1/2 up to c = 0. Multiplying 
the result by a factor of two gives the desired result for the 05^(212) case. 

We now get back to the untwisted sector with Cgg- = 2, and discuss the finite size spectrum. 
Excitations are obtained by changing the number of roots or shifting the roots. Parametrize the 
numbers of roots as 

L , , 

n,/7 = ^ - (90) 

and suppose the Bethe integers are submitted to a global shift given by miji. One then finds 
(similar formulas appear in |9j) 



. _ C\ f ni + rajj Y 1 / mi + m jj ^ ^ 

+ i(Q)%^(!=i^)' (91) 

For reference it is useful to recall the finite size spectrum of the XXZ antiferromagnetic chain 
with A = — cos ^0 '■ 

io\ 1 m? 



1 - - — + — - (92) 

vr / 2 1 _ & 2 ^ ' 

We observe first that excitations with n/ = nn = n and m/ = mn = m in the OSp chain are 
equal to twice the excitations in the XXZ chain, a result that can be checked directly at the level 
of the Bethe equations. 

The OSp quantum numbers are given by 

^ _ ni - mi 
2 

= (93) 
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so the sub-component of the spectrum arising from excitations mj = mjj = reads 



X(n„n„)(0,0) = - J j (^') + (94) 

We see on (I9ip that the finite size spectrum is determined by two bosonic fields with related 
coupling constants. When ^ — t- 0, one of the coupling constants vanishes, and the theory reduces 
to a non compact boson and a compact one at the Dirac radius. 

We now consider, like in the case ^ = 0, a staggering of the OSp{2\2)g chain (adding ±iA 
to U,T in the Bethe equations), in order to induce a massive deformation. We now use Fourier 
transforms defined as 

m = [ ^e^^^/2^/(C/) (95) 



The equivalent of equations (|M|) , is now 

I, cosA;A/2^ ^ h ^ h 

+ = + + + • • • ^''^ 

where we have not written the interaction with the other strings, and the kernels are given by 

coshfc(f - 2) - 1 



= ^ii,ii 



2sinh/c(| — l)sinhA; 



coshfc(| — 1) — cosh A; 
= = 2sinhMf-l)smhA: ^^^^ 

The dominant pole in the source term is a A; = zvr, so the mass term goes as M oc e"^'^/^^, and 
the free energy will be a function of the product M/T as usual. We have then to match analytical 
properties of the free energy deduced from the lattice equations with the relation between the 
physical mass and the bare mass (related with the staggering). A few trials seem to leave only 
one option: a perturbation of the type 

A = j fx^ + Acos/3i(/)ie^^2<^2 + ^g-'^iM^ (98) 



=1,2 



with the constraint 



/3? + /3| = 47r (99) 

Eq. (I99p follows from the finite size lattice spectrum, where pure 4>i excitations are obtained with 
mi = rrijj = and n/ = ±njj. To proceed, we observe that the action in ([98]) leads to a free 
energy expanding in even powers of AV^^'^^'/" = AV^^'/" (where we used (|99l)). It follows 
that 

Moc(AV)'^ (100) 

Meanwhile, the perturbation induced by the staggering on the microscopic hamiltonian is pro- 
portional to e~^, so we get the basic identification ^ = or 

^ = L 

Svr 27r 

^ = i-A (101) 

Svr 2 27r ^ ' 
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This is in the finite size spectrum indeed. Pure (pi excitations correspond to n/ = —nu, pure 
4)2 excitations to nj = n//. When — )• 0, we get in the action, to leading order, 6*^^"^^ and 
^-2i\/An(f)2 ^ that is, an a'^^ action, which is known [1] to be equivalent to the pure fermion terms 
in (I73p after twisting (to et c^gf = 2). The (pi field, meanwhile, must correspond to the non 
compact <j) field in ()73p . We recover the expectation that the non compact direction corresponds 
to excitations nj ^ nu then. 

We note now that the model (|98|) is a complex version of the model dubbed c[^^ in [38], and 
well known to be integrable. Our analysis gives as a by-product the free energy of this theory: 
similar arguments as in earlier sections about the sector symmetric under I -H- II exchange shows 
that 

^(C« - complex) = 2^(sG,^ = f] (102) 



L' ' ' ' L \ ' Sir 47r, 

where (1 1011) must be used to obtain We do not know how such a relation could be proven 
directly. 

Another intriguing observation is that the kernels for the I — I and I — II scattering, which 
read (compare with the isotropic case (|67p recovered by letting ^ — t- oo and identifying k = x) 

Id \ r coshff - 2)k - 1 

ide vri 2sinh(|-l)A:sinhA; 

Id If cosh(f — 1)A; — coshfc , 

-im = J "'^idif^W^' 

(here 9 is the physical rapidity, 9 = -^U, obey the relation 

Id If sinh(;^ — l)k , 

id9 '''' ttJ 2sinh(^ - i)fccosh| ^ ^ 



This is exactly the sine- Gordon kernel at 



1 - - (105) 



A more detailed analysis shows that combining holes of type I and type // gives excitations 
whose scattering is exactly given by the S matrix for the q-deformed GN theory in fl2] . It is also 
possible to use the limit g — )• 1 to try to make sense of the diverging kernels encountered in the 
isotropic case. This will be discussed elsewhere. 

Note finally that when 7 = f , the theory seems entirely free: no interactions remain in the 
Bethe ansatz, and Zjj = Zjjj = 1. This corresponds meanwhile to the case (3i = (^2 = V^tt, 
which will also be discussed elsewhere. 



8 Application: flow in the Random Bond Ising model 

In a very interesting paper, Mussardo et al. [39] argued that the flow into the random bond 
Ising model could be described by an 0{N) massless scattering in the 'limit' — ?• 0. The 
S matrix proposed in [39] involves left movers in a copy of the fundamental representation of 
0(A), right movers in another copy, with Srr = Sll = Srl- This is the obvious structure one 
might guess from experience in the flow from 0(3) to SU{2)i WZW model. Such an S matrix 



28 



Figure 13: TBA for the susy version of the massless flow into the random bond Ising model. 




Figure 14: TBA for the susy version of the massless flow into the random bond Ising model: an 
equivalent form after folding. 

structure has also been considered by P. Fendley in [lO]. There, the question being studied was 
the existence of a massless flow into the 0(2A'')i WZW model. The S matrix proposed had a 
similar structure, Srr = Sll = Sgn (the case of general N is considered, with particles in the 
fundamental representation). The S matrix for Sir in [50] required some additional CDD factors 
to get rid of spurious poles in the LR channel. It was then found to describe a flow from the 
0{2N)/0{N) X 0{N) Sigma model at = vr into the 0{2N)i WZW model. 

Of course the regime studied in [40J concerns N > 2, and has a different behavior from the 
one in [39] since relevant and irrelevant directions are switched as crosses two. We thus cannot 
use the results in [30] analytically continued in A^ to see what happens in the random bond Ising 
model. 

The natural strategy to make progress is then to turn to a supersymmetric version |i4lj of 
the argument in [39j, replacing the A^ — )• limit by OSp{2\2). The goal is then to find a RG 
trajectory ending up at the 05^(212)1 WZW model in the IR, and originating form a non trivial 
UV fixed point. The proposal in [39] immediately extends to Srr = Sll = Srl given by the S 
matrix in section 4, with R particles having dispersion relation e = p = ^ e^^, and L particles 

Note that, in contrast with the case discussed in [30], no CDD factor is a priori necessary, 
since there are no poles in the GN S matrix any longer. 

The TBA in the Neveu Schwartz sector is thus obtained by an immediate generalization of 
the results in the previous section, where instead of having a massive node on the left, we have 
a L and a R moving node, identically coupled to the body of the diagram, with no coupling 
between them. See figure I13I14I and the appendices. It is clear that this TBA converges to the 
correct theory in the IR. The question is, where does it originate from? 

Let us write the TBA explicitely. We give to the two grey nodes labels 0,0, and to the other 
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nodes labels 1,2, We have thus 



eQ(^) = ^e-'-2Ts^ln{l + e~^^/^ 
ei = —Ts ★ 



€2 = -Ts -k 



In (l + e-"o/^) +ln(l + e~'o/T^ +ln(l + e'^^/^^ 
In fl + e-'^/^] + In fl + e~''^/^] + ... 



and 



(106) 

J In (l + e-^«/^) + ^e"^ In (l + e-^o/^) (107) 



where s = t, r-s and * denotes convolution. 

27r cosh 

This is quite similar to the TBA for the 0(3) sigma model at = vr; in fact, the only difference 
is some factors of two in the coupling to the source nodes. 

The game now consists in solving these equations numerically (after some truncation) and 
studying the IR limit. This is a long technical story delineated in the appendix. The main 
conclusion is that the TBA develops a singularity at some value of M/T, and that this value is 
very large, and probably infinite. Such intriguing behavior was studied in |42j for instance, where 
it was interpreted as a sign that the theory one is trying to define encounters some sort of phase 
transition. For the examples discussed in [32], the singularity remained at a finite value of M/T, 
so the vicinity of the IR fixed point at least was well defined. In our case, it seems the singularity 
might occur right near the IR fixed point, so the WZW theory sits on the border of a domain of 
instability in the massless region. Meanwhile, we have found that adding a CDD factor in the 
LR scattering seems to give a meaningful TBA, provided one chooses 

/ = tanh2(0/2 + i7r/4) (108) 

Note that this factor has a double pole at 9 = i/2, which leads to a 'monstron' with complex 
mass 

Mm = Me*^/^ (109) 

similar to the one encountered in the study of the massless flow from dilute to dense polymers 
\A'6\ 157] . The TBA now gives Cgg = 1 in the UV, but it is not clear whether one should add to 
it a monstron contribution [H] , [37] . We hope to present a more detailed study of this question 
elsewhere. 

9 Conclusion 

On the positive side, we have, in this paper: 

• built the basic TBA for OSp{2\2) symmetric integrable models 

• found evidence that the S matrix of [12] for the GN model is the correct one 



uncovered an intriguing duality between different quantum fleld theories, somewhat gener- 
alizing to the massive case the results from [5l [6] 
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• uncovered yet another duality between the (q-deformation of) the 05^(212) sigma model 
and some integrable two bosons theories 

• shown that the proposal of [39J for the flow into the random bond Ising model is not correct 
We have also uncovered some troubling features that seem very generic. In particular 

• the S matrix of the OSp{2\2) sigma model appear to be a very intriguing object, where 
some of the usual principles might not apply. It is not known up to now, even though the 
TBA has been extracted from a lattice regularization. It is also not known what kind of 
'sausage deformation' of the sigma model corresponds to the q-deformation of the lattice 
equations. 

It also remains to see how our analysis might fit into the more general scheme recently 
developed in the AdS/CFT enterprise. We note in this respect that the problems discussed in 
this paper could a priori be tackled as well with another form of the Bethe ansatz, corresponding 
to choosing a different Dynkin diagram. This question was briefly alluded to in [32j. We leave 
this question to further study, and simply notice that the Y system for the GN model reads 
(setting as usual Y = e~^^'^) 

Yo{e-m/2)Yo{e + iTT/2) = {l + Yi{e)f 
Y^{e-m/2)Y-^{e + iTi/2) = {l + Yi{e)f 
Yi{e-m/2)Yi{e + m/2) = (I + Yo{e)){l + Y^){e)){l + Y2{e)) 

Yi{e-m/2)Y{e + iTi/2) = {i + Y-i{e)){i + Y+i){e)) (iio) 

This differs from the one for the 0(3) model only by a power of two appearing in the first two 
equations. It would be most interesting to match this Y system with some of the general results 
uncovered recently |l5l HH |3Z] . 

We would like to end this paper with a speculation. We have identified the q-deformed 
OSp{2\2) sigma model with a complex version of the theory discussed in [38]. Now, for the 

real case, V. Fateev conjectures that the theory is dual to the complex sinh-Gordon model 
(CSG). Extending (maybe too naively) his conjecture to the complex case, we find that our q- 
deformed OSp{2\2) sigma model should be equivalent to the complex sine-Gordon model with 
action 

1 [ d^xdtiX 2u,i2 



A=- — m-nlyl + counterterms (HI) 

2 J I - ^|y|2 

In particular, we obtain that the pure 05^(212) case should correspond to 

A = - — - — p— ^ — mnlxP + counterterms (112) 

2y i-TT\xr 

Of course, the coupling is now so big in the CSG that it is not clear what the action means. Yet, 
another important observation in favor of this proposal concerns conserved quantities. It is well 
known that the CSG model admits independent conserved quantities of arbitrary spin [48J . On 
the other hand, the S matrix encodes conserved quantities for instance, expanding the sine- 
Gordon S matrix (|104p . the poles at k = (2n + l)z7r correspond to the local conserved quantities, 
which act as e(2"+i)^ on particles, and thus are present only for odd spin. Meanwhile, in the 
conjectured S matrix (I103p . there are poles at all k = nvr, leading also to conserved quantities 
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acting as e " . These poles disappear in the product Zj^Ziji but are otherwise present, as would 
be required for the CSG. 

Interestingly, the question of S matrices for the complex sine-Gordon model is far from being 
settled [50l |51] . Our results provide an intriguing light on this question, indicating a spectrum 
quite different from the semi classical proposals. We will discuss this issue in more details else- 
where. 
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Appendix A: Numerical study of the massless TBA 
A.l 0(3) sigma-model with ^ = tt: a reminder 

We first consider this case as a warmup. The infinite set of TBA equations is well known (see 
also Fig. USD: 

eo = re72-s*ln(l-he"'^) 

eo = re-V2-s*ln(l + e-^0 (Ai) 
€i = -s ★ In (1 + e~'o) - s * In (l e"'5) _ g * In (l e"'^) 

= -s*ln (1 -h e"''-i) - s*ln (1 e"''+i) for « = 2,3, ... 

with s = cosh 6 ' have rescaled the e's with respect to the main text, and set r = M/T. 
The dimensionless free energy per unit length F[r) is calculated as 

In the UV and IR limit it is related to the central charge as F = — c/12. 

••• -o o 



Figure 15: TBA for the 0(3) sigma-model with 9 = n 

The above set of equations has to be truncated to be accesible to numerical investigations. 
Let us denote the total number of nodes with N. The expected values of the central charge are 

_ 6 _ 6 

''^^ " iVT2 - 1 - (^^2)(A + 3) 

The function F(r) is plotted for two different values of N in Fig. [THl 

Convergence as a function of N is clearly rather quick. To see this better, consider F as a 
function of N for different fixed values of r. The results for r = 0.1, 1, 10 are shown in Fig. [T71 
One can see that there is a well-defined A — t- oo limit, the leading correction is 0(1/A). 
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logr 

(a) iV = 6 (b) AT = 10 

Figure 16: F(r) in the truncated versions of the cr-model TBA 
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Figure 17: F as a function of N for r = 0.1, 1, 10 in the 0(3) sigma-model with 9 = ir 
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Figure 18: TBA for the susy version of the massless flow into the random bond Ising modeh an 
equivalent form after folding. 

A. 2 The conjectured massless flow in the 05'P(2|2) case 
The infinite set of TBA equations given in the text becomes: 

eo = reV2-2s*ln(l + e~'i) 
eo = re-V2-2s*ln(l + e-^i) 

= -s * In (1 + e"'») - s * In (l + e"'o) - s ★ In (l + e"'^) 
ei = -s*ln(l + e-'»-i) -s*ln(l + e"''+i) for i = 2,3,... 

with s = cosh d ' "^^^ Only difference compared to the case of the o"- model is a factor of 2 in the 
equations for eo and eg. However, this factor dramatically changes the behaviour of the TBA. 
Similar to the cr-model, we first study the truncated systems with nodes. 

The unfolded version of the TBA is represented by the graph in fig. [TSl which only contains 
simple links. Given this structure of the TBA, only those diagrams can describe a meaningful 
flow for which the largest eigenvalue Xmax of the incidence matrix is less than 2 [51]. If this 
condition is satisfied, in the UV limit one observes the usual pattern of platens in the functions 
Li = log(l + e"*^') in the central region \9\ < logr. On the other hand, if X^ax > 2, then one 
expects divergences in the UV limit This is the case for the diagram [THl even for finite A^. 
Nevertheless the calculation for the central charge is well-defined in the IR and leads to 

^ 3N-1 

We performed numerical calculations for different values of A^, the results for A^ = 6 and 
A^ = 10 are shown in Fig. [191 One observes that the truncated equations work fine in the 
IR and produce the predicted central charges. However, the iterations become unstable (the 
pseudoenergies and also F diverge to minus infinity) whenever r is below some critical (A^- 
dependent) value. One finds moreover that rcrit increases with A^. 

We determined rcrit for different values of A^, the results are shown in Table 1. and Fig. [501 
The numerical iterations become more time demanding as A^ is increased, therefore we had to 
content ourselves with an accuracy of Ar = 0.1 for the larger values of A^. 

At this point it is not clear what happens with Vcrit as N is increased further. There are two 
possibilities: 

1. There exists a limiting value for rcrit as A^ — t- oo. In this case the infinite set of TBA 
equations would be meaningful in the interval r = rcrit ■■■ oo. However, one would see 
unusually high effective central charges, since already for A^ = 20 and r = 4 one has 
Ceff ~ 10. 

^Xmax may be equal to 2, but then this should hold for every finite N. In this case the pseudoenergies may be 
infinite in the UV limit, whereas the elTective central charge stays finite. 
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Figure 19: F{r) in the truncated version of the massless TBA ()A2p . The dashed hne corresponds 
to the analytical prediction in the IR. The filled circles show the critical values of r where the 
TBA equations become unstable. 



2. Vcrit — )• oo as — )• oo. 

To gain more insight, one can consider F as a function of N for different fixed values of r; 
the results are shown in Fig. [2TJ For r < 5 we observed, that the iterations become unstable 
whenever N reaches a critical value Nc < 100. These points are marked with solid squares in 
the figure. For r > 5 we did not find the critical values of A^, because we had to restrict the 
calculations to A < 100. However, from the figure it is clear that even for r > 5 there is no 
convergence. In fact one observes a linear gorwth in N for every value of r, suggesting that even 
for r > 5 there always exists a (possibly quite large) critical N. This means in turn, that in the 
A — 7- oo limit one has rcrit oo. 

We tried the improve the TBA by introducing alternative truncation schemes for the infinite 
diagram, for example putting a "fork" at the end. However, none of these tricks works to stabilize 
the iterations for small r, because Xmax will still be greater than 2. 

Based on the evidence presented above, we conclude that the infinite set of equations (|A2p 
cannot be solved for any finite r. 



Appendix B: Modifications of the massless TBA 
B.l CDD-factors 

It is a natural idea to add CDD factors to the conjectured S-matrix and to consider the resulting 
TBA equations. CDD factors in the LL and RR channels would alter the IR central charge, and 
therefore we can exclude them. The only possibility is to introduce CDD factors for Srl and 
Slr- 

A simple factor of 

^ sinh 6 — i sin an 
sinh 6 + i sin an 

introduces a coupling between the massless nodes in the TBA equations 



eo = - 2s ★ In (l + e'^^) - (^^ * In (l + e'^o) 

eg = re~V2 - 2s * In (l + e~'') - ★ In (l + e~'«) 

ei = -s * In (1 + e~'°) - s * In (l + e~'5) - s ★ In (l + e^'^) 

ei -s*ln(l + e-'"-i) -s*ln(l + e~''+i) for i = 2,3, ... 



(Bl) 
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where 

_ 1 2 cosh 6 sin air 
2tt sinh 9"^ + sin^ air 

The function ipa has the same sign as s = cLh e if « > 0, therefore it makes the situation 
even worse. The only possibiUty to "cure" the TBA is to choose a < 0. In this case there is 
a pole in the physical sheet at ^ = ia, which should be explained. However, the TBA makes 
sense perfectly, irrespective of the value of a. In fact, a can be chosen arbitrarily, because it only 
determines the crossover scale M. Choosing a = 1/2 one obtains 

/ = tanh^(6'/2 + i7r/4) with (p 



IT cosh 9 



The folded version of the TBA is then represented by the graph below, where the vertical line 
with the arrows represents ,, negative" coupling. 




The UV central charge cannot be obtained in the usual way, because there are still no finite 
platens in L = log(l + e~^) and in the central region one has = — oo for every node. Instead 
one can consider the truncated diagrams with N nodes and let A/" — >■ oo. This way one obtains 

Ceff = l (UV) Cgff = 3 (IR) 

The above results were obtained by numerically solving the Y-system equations 

Yo{9 - m/2)Yo{e + m/2) = (1 + Yi{e)f ^ 



+ 

Y^{9 - m/2)Y^{9 + m/2) = (1 + ^22) 

Yi{9 - m/2)Yi{9 + m/2) = (1 + Yi,{9)){l + Y^{9)){1 + Y2{9)) 
Yi{9 - m/2)Yi{9 + m/2) = (1 + y^-i(^))(l + 

with constant numbers yi and evaluating the central charge with the usual rules. We do not know 
of an exact analytical proof. However, we can give an approximate solution, as follows. 

Let us consider a truncated diagram with N nodes denoted by 0, 0, 1, . . . {N — 2). In the UV 
limit the variables yi = e^*^' take the values 

^ / 27(Ar-2)2 ^9(iV-2)2/ /27r(7V-2-j)\ \ 

in the central region. On the other hand, for — > oo the constant values are approximately 

4(A^-2)2 2(Ar-2)2/ /7r(Ar-2-z)\\ 

2/0 = y-o^ \2 ' Vi^ ^2 [l-cosi^ ' jj (B4) 

Actually, in this latter case the exact solution is known, because this corresponds to the Y-system 
equations of the An diagram with n = 2N — 3. However, the exact solution is not needed. The 
effective central charge only depends on the way the variables yi approach infinity for N oo 
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in ()B3jl and ()B4p . Substituting the above values into eq. (5.3) of [53] yields c = 1 in the N ^ oo 
limit. 

We performed numerical simulations of the truncated equations (IBlh with different values of 
N between = 10 and = 40 and extrapolated the results to — )■ oo applying a polynomial 
fit in The results are shown in Fig. [22l 

We note that a simple factor of tanh(0/2 + z7r/4) is not sufficient to stabilize the TEA, because 
in this case one still has A > 2 even for the truncated diagrams. 

B.2 Slr being a pure CDD factor 

In this case (which was in fact discussed earlier in the literature [SU]) the theory posseses two 
independent OSP(2\2) symmetries corresponding to LL and RR scattering. The TEA we are 
dealing with is now represented by the diagram 

... -o o o o o- ... 

(-) 

... -o o o o o- ••• 

The vertical line connecting the two inifinite chains corresponds to the CDD factor. A simple 
factor of 

Slr = t&nh(e/2-i7r/ 4) 

(which proved to be the correct choice in [55] and |56j ) does not work in this case, because once 
again the TEA will be oversaturated. The simplest possibility is to choose the inverse of the 
latter: 

Slr = t8inli{e/2 + in/ 4) 
In this case the dashed link corresponds to 

1 



2-K cosh 6 



The TEA now gives Cg// = in the UV (the truncated diagrams have finite positive Ce//, which 
tends to zero if A^ — )• oo), which is a probably meaningless result. 
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N 




4 


0.32175 


6 


1.305 


8 


2.016 


12 


2.95 


20 


4.1 


30 


4.7 


50 


5.5 



Table 1: rcrit for some values of N 
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Figure 20: Critical value of r as a function of N 
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Figure 21: F as a function of N for different fixed values of r. The solid square represent the 
critical values of N, where the iterations become unstable. 
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Figure 22: The evolution of the effective central charge in the massless TBA with a CDD factor. 
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